Probing single charge fluctuations in a semiconductor with laser spectroscopy on a quantum dot 
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We probe local charge fluctuations in a semiconductor via laser spectroscopy on a nearby self-assembled 
quantum dot. We demonstrate that the quantum dot is sensitive to changes in the local environment at the single 
charge level. By controlling the charge state of localized defects, we are able to infer the distance of the defects 
from the quantum dot with ±5 nm resolution. The results identify and quantify the main source of charge noise 
in the commonly-used optical field-effect devices. Based on this understanding we achieve routinely close-to- 
transform-limited quantum dot optical linewidths. 

PACS numbers: 73.21.La and 78.67.Hc 



Condensed matter systems, notably quantum dots in III-V 
semiconductors and color centers in diamond, are very attrac- 
tive as the building blocks for quantum light sources |Q] and 
spin qubits J2j. For instance, an InGaAs quantum dot is a 
robust, high repetition rate, narrow linewidth source of on- 
demand single photons and polarization-entangled photons, 
properties not shared by any other emitter. In the future, the 
demands placed on the quality of the single photons will in- 
crease. For instance, the creation of remote entanglement via 
photon interference and associated applications as a quantum 
repeater require Fourier-transform-limited single photons, i.e. 
wavepackets with a spectral bandwidth determined only by 
the radiative lifetime. This is hard to achieve in a semicon- 
ductor. On the one hand, a quantum dot is extremely sensitive 
to the local electric field via the Stark effect [3 4] leading to 
a stringent limit on the acceptable charge noise. Charge noise 
can also lead to spin dephasing J5j|6). On the other hand, 
phonons in the host semiconductor can lead to dephasing JTJ. 
However, at low temperature and with weak optical excita- 
tion, phonon scattering is suppressed in a quantum dot by the 
strong quantum confinement J8j|9], and the remaining broad- 
ening arises from relatively slow fluctuations of the environ- 
ment leading to spectral fluctuations ifTUl . Transform-limited 
lines have not been routinely achieved, with typical optical 
linewidths a factor of at least 2 or 3 above the theoretical 
limit iTTOlfPJl . While spectral fluctuations in self-assembled 
quantum dots have been investigated with non-resonant exci- 
tation lfT4l [131 . their origin in the case of true resonant ex- 
citation is not known with any precision and are potentially 
complex with contributions from a number of sources. Fur- 



thermore, spectral fluctuations are a common feature in con- 
densed matter systems, arising also in diamond [16], semicon- 
ductor nanocrystals [17] and nanowires [18|. 

We report new insights into local charge fluctuations in a 
semiconductor. High resolution laser spectroscopy on a sin- 
gle quantum dot is used as an ultra-sensitive sensor of the local 
environment. We observe single charge fluctuations in the oc- 
cupation of a small number of defects located within ~ 100 
nm of the quantum dot. We control the occupation of these 
close-by defects with an additional non-resonant excitation. 
Once the defects are fully occupied, there is a strong sup- 
pression of the charge noise. This understanding is tested in 
a new heterostructure in which the fluctuators are positioned 
further away from the quantum dot. As predicted by our mod- 
elq, this change reduces significantly the quantum dot opti- 
cal linewidth, making the observation of close-to- transform- 
limited linewidths routine. 

The InGaAs quantum dots are embedded in a Schottky 
diode fT9ll20"l , Fig. 1(a). They are separated from an n + back 
contact by a dtun = 25 nm thick GaAs tunnel barrier. Di- 
rectly on top of the dots is a capping layer of thickness d cap , 
30 nm in samples A and B, followed by a blocking barrier, an 
AlAs/GaAs superlattice: c?sl = 120 nm in sample A, 240 nm 
sample B. Sample C has d cap = 150 nm and rfgL = 240 nm. 
Samples B and C were grown under identical conditions. The 
samples are processed with Ohmic contacts to the back con- 
tact, grounded in the experiment, and with a semi-transparent 
gate electrode on the surface to which a gate voltage, V 9 , is ap- 
plied. Laser spectroscopy is carried out on the charged exciton 
X x ~ at 4.2 K by focusing the linearly -polarized output of a 1 
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FIG. 1. (a) Band diagram of the devices, (b) The optical set-up for 
AR/R measurements, (c) AR/R versus gate voltage for constant 
resonant laser wavelength (951.1 150 nm) and power (1.0 nW) for a 
quantum dot in sample A (d cap = 30 nm) both without (black) and 
with (red) P — 325 nW of 830 nm laser light. The inset shows the 
resonance position versus V g . The Stark shift depends linearly on 
voltage away from the plateau edges; the Stark parameter decreases 
by only 10% at P = 325 nW. (d) Color-scale plot (linear scale, blue: 
0.061%; red: 0.61%) of AR/R versus non-resonant laser power P. 



MHz linewidth laser to a ~ 0.5 /im spot on the sample surface. 
The power of the resonant laser is ~ 1 nW to avoid power 
broadening. The key advance here is to illuminate the sample 
simultaneously with a weak non-resonant source at 830 nm, 
Fig. 1(b), with power P. Resonant excitation of the quantum 
dot is detected either with differential reflectivity AR/R ET1 
including a filter to reject the 830 nm light Fig. 1(b) or with 
resonance fluorescence exploiting a dark field technique |fT3l . 
The integration time per point is typically 500 (100) ms in 
AR/R (resonance fluorescence). Spectra are recorded either 
by sweeping V g (changing the detuning via the Stark effect) 
or by tuning the laser. 

Fig. 1(c) shows typical laser spectroscopy results both with- 
out and with "high" non-resonant excitation, P = 325 nW. 
In both cases, the absorption lines are close to Lorentzians 
with linewidth 2.5 /xeV. The radiative lifetime at this wave- 
length is 800 ps [22 1, implying transform-limited linewidths 
of 0.8 /xeV, a factor of 3 smaller than observed in the experi- 
ment. Other groups achieve similar linewidths lfToUT3ll . The 
main effect of the non-resonant excitation is to shift the reso- 
nance to more negative voltages, in this case by AV g = —80 
mV, for the same laser wavelength, equivalently a blue-shift 
of AE — 60 //eV for the same gate voltage, Fig. 1(c). Fig. 
1(d) shows AR/R over 4 decades of P. Remarkably, the dot 
evolves from the low-P region (single Lorentzian line inde- 
pendent of P) to the high-P region (single Lorentzian line 
shifting monotonically with P) via a series of steps. These 



steps occur rather abruptly, over just a decade in P. For this 
particular quantum dot, 4 steps (equivalently 5 AR/R lines) 
are observed. The energy separation of the lines varies from 
about 4 to 10 /ieV, and the linecuts, Fig. 2(c)-(e), show that 
within each line there is also a sub-structure. The observa- 
tion of these absorption steps and their behavior as a function 
of the control parameter P constitute our main experimental 
discovery. 

We find that the P = and P = 100 nW behavior are 
very similar for all dots. Also, the intrinsic properties (radia- 
tive lifetime, Stark shift, Coulomb shifts on charging) are all 
broadly similar. Despite this, the transition region is highly 
dot dependent. The number of steps lies typically between 3 
and 6; the energy separations between the lines lie between 
~ 4 and 20 /xeV (for sample A) with each quantum dot hav- 
ing its own unique "finger print" in the P-dependence. We 
therefore look for an explanation of the steps in terms of the 
dots' environment, i.e. a mesoscopic effect. 

Our hypothesis is that nonresonant excitation creates holes 
at the capping layer/blocking barrier interface, Fig. 1(a). 830 
nm light creates electron-hole pairs in the wetting layer. The 
electrons relax rapidly to the back contact, the holes to the 
capping layer/blocking barrier interface where at low temper- 
ature they can be trapped, creating a positive space charge in 
the device. The trapped holes mean that the same electric field 
is achieved at the location of the quantum dot only by apply- 
ing a more negative voltage to the gate, consistent with Fig. 
1(c). At large P, a 2D hole gas is formed, and the shift in 
V g of the optical resonance allows the hole density Nh to be 
estimated. For intermediate P where we observe the steps, 
the hole density can be estimated for sample B to be ~ 10 10 
cm~ 2 , similar to reported values at the metal-insulator transi- 
tion 11231 . The steps arise in the localization regime. In partic- 
ular, the steps reflect a change of just one hole in occupation of 
the localization centers close to the dot. Quantitatively, occu- 
pying a localization center immediately above a quantum dot 
at d cap — 30 nm changes the electric field by —1.50 kVcm -1 
(taking into account the image charge in the back contact), 
shifting the optical resonance by 20 /ieV via the Stark shift. 
This corresponds closely to the maximum observed step sepa- 
ration. This, and the agreement with our simulations (below), 
justifies our hypothesis. Smaller steps arise from the occupa- 
tion of localization centers which are laterally displaced. 

Our interpretation leads to two immediate results. First, 
the location of the energy line of the quantum dot is a direct 
measure of the number of charges stored directly above the 
quantum dot. In the low-P regime, the quantum dot senses 
the nearby environment with single charge resolution. Sec- 
ondly, the number of steps observed equals the number of 
holes which can be trapped above the dot, 4 for the dot in 
Fig. 1(d). 

We underpin our experimental results with a Monte Carlo 
simulation of the effects of occupying an array of valence 
band localization centers at the capping layer/blocking bar- 
rier interface. We take an array of localization centers all dis- 
tance d cap above the quantum dots but at different locations 
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FIG. 2. (a) AR/R versus P for a quantum dot in sample A, as in 
Fig. 1(d). (b) Monte Carlo simulation with 4 hole localization centers 
located above the dot with n = (32.0, 15.4, 15.7, 48.8) nm, a t = 
(5.0, 1.5, 1.5, 1.0), T = 2.5 /ueV, T L = 1.0 /im and N = 2, 500 
(parameters described in the text), (c)-(e) Line cuts at four different 
values of P showing experimental data (black) and simulation (red), 
(f) Lateral location of localization centers with dot at r = 0. 



Vi = (ri,9i) within the 2D plane. We position by hand a 
small number of localization centers, between 1 and 4, each 
with r < 50 nm. Additionally, we take a full 2D array of ran- 
domly placed defects with 2D density A^zj- The occupation 
of a defect changes the local electric field at the quantum dot 
and hence the absorption spectrum via the Stark effect. This 
is calculated by, first, calculating the additional electrostatic 
potential; second, the associated electric field; and third, the 
energy shift of the exciton via the Stark effect. The Stark shift 
from the vertical electric field is calculated from the measured 
Stark effect, i.e. from the Vg-dependence of each particular 
quantum dot (modeled as a permanent dipole moment in an 
electric field ll24l ). The lateral electric field component can- 
not be accessed directly in the experiment but the effects are 
smaller: we assume that there is no linear term (i.e. no perma- 
nent dipole moment in the lateral plane) and that the quadratic 
component scales with the fourth power of the wave function 
extent of the quantum dot ground state which is known rea- 
sonably well lt24l . The localization centers i are each occu- 
pied with a probability a.ip which rises with p, the control 
parameter in the simulation, until it reaches 100%. cti can 
vary from center to center and expresses the relative probabil- 
ity of occupying a particular center. With a full 2D array, on 



depends on n through a Gaussian function with full-width-at- 
half-maximum (FWHM) Tl which describes the spatial ex- 
tent of the non-resonant beam focus. For the defects directly 
above the quantum dot, the a, are treated as fit parameters. 
For a fixed defect distribution and for a given p, we occupy 
the defects with a random number generator; from this charge 
distribution we calculate the net Stark shift, and at this energy 
we place a Lorentzian absorption spectrum with FWHM T. 
The process is repeated A" times, keeping the defect distribu- 
tion constant but each time creating a new charge distribution 
with the random number generator. The whole procedure is 
then repeated as a function of p. We model the experiment by 
relating p linearly to the control parameter P. 

Our simulation reproduces the steps in the absorption spec- 
tra as a function of P for sample A, adding considerable 
weight to our assertion that the charge fluctuations arise from 
trapped holes at the capping layer/blocking barrier interface. 
The exact energy steps turn out to be very sensitive to the loca- 
tions Ti of the localization centers. (The dependence on 0i is 
much weaker.) We can match the energies of the steps, their 
P-dependence and the substructure within each step with a 
set of Ti and N 2 d = 0. However, we need to depart from 
ci, = 1 to reproduce the relative intensities of the various 
lines. Fig. 2 shows the result of this procedure: the Monte 
Carlo simulation, Fig. 2(b), reproduces the main experimen- 
tal features, Fig. 2(a). Furthermore, the line-cuts at specific P 
are in very close agreement with the complicated experimental 
spectra, Fig. 2(c)-(e). The localization center array is shown 
in Fig. 2(f). The different a, presumably reflect some con- 
nectivity between the localization centers such that a "deep" 
one is much more likely to be occupied than a "shallow" one. 
The energy shifts on adding holes to these defects one by one 
are so sensitive to the set of Vi that the random error on each 
ri is as small as ±5 nm. In this sense, the experiment pro- 
vides ~ A/100 spatial resolution in the spacings between the 
localization centers. 

We have attempted to reproduce results such as those in 
Fig. 2(a) with just a random distribution of localization cen- 
ters, N 2 d 7^ 0- The large net shift between P "low" and 
"high" pins down N 2 d to ~ 10 10 cm -2 . For this N 2 d, the 
Monte Carlo simulations predict only in very rare cases 3 — 5 
steps yet this is the standard experimental result. Furthermore, 
in the simulation, each line has a strong P-dependence, not a 
feature in the experiment. In the simulations, the only con- 
figurations which describe sample A are those with a cluster 
of localization centers immediately above the dots with oth- 
erwise a sparse distribution for r < 100 nm, an extremely 
unlikely outcome with a random distribution of localization 
centers. The conclusion is that the localization centers are not 
randomly distributed in the 2D plane. Instead, the dot itself in- 
duces the formation of a small number of localization centers 
directly above it. The mechanism for this is likely to be the 
strain field which extends beyond the quantum dot in com- 
bination with roughness at the capping layer/blocking layer 
interface. 

Sample A has T = 2.5 /ieV, well above the transform- 
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FIG. 3. Resonance fluorescence from a single quantum dot in sample 
B (1.0 nW at A = 962.2500 nm) at P = with integration time 
per point 0.1 s in (a), 50.0 s in (c). Monte Carlo simulation with 
N 2D = 1.0 x 10 10 cm -2 , T = 0.8 ^eV, F L = 10.0 /im, p = 4.4% 
(to represent P = 0) with N = 10 in (b), N = 5, 000 in (d). 



limit. As described above, this is unlikely to be related to 
fluctuations in a 2D array of localization centers at the capping 
layer/blocking barrier interface. The origin of this broadening 
is not known precisely but there are hints that it is related to 
the surface of the device. We switch to sample B which clearly 
demonstrates the consequences of a fluctuating 2D array. Fig. 
3(a) shows resonance fluorescence from a single dot in sam- 
ple B. At P = 0, the FWHM is comparable to those of dots in 
sample A but there are large fluctuations in the signal which 
are not reproducible from one spectrum to the next. The fluc- 
tuations disappear only when we integrate for more than 50 s 
per point, Fig. 3(c), demonstrating that they have a component 
at very low frequency (sub-Hz). A characteristic feature is the 
rather abrupt turn on at negative detunings and the abrupt turn 
off at positive detunings. Turning on the non-resonant exci- 
tation reveals also a series of steps, as in Fig. 1(d), and at 
"high" P, this sub-Hz frequency component is eliminated. We 
interpret the P = results with the Monte Carlo simulations, 
Fig. 3(b), with the hypothesis that the /ieV-scale fluctuations 
in Fig. 3(a) arise from fluctuations amongst a large number of 
localization centers all with r < 100 nm. With this hypothe- 
sis we can reproduce the experiment, Fig. 3(b),(d), provided p 
is small, i.e. the defects are each occupied with small proba- 
bility. The defects (two in this case) directly above the quan- 
tum dot are therefore unlikely to be occupied. Only a small 
fraction of the available configurations are occupied within 
the measurement time, leading to the changes in spectrum to 
spectrum. Fig. 3(b) reproduces the abrupt turn on/turn off of 
the spectrum, the FWHM, and the characteristic energy split- 
ting between the sub-peaks using N2D = 1-0 x 10 10 cm -2 , 
p = 4.4% and N = 10. Significantly, the jagged nature of the 
spectra in Fig. 3(a) can only be reproduced with a small ho- 
mogeneous broadening, V = 0.8 fisV. This is evidence that on 
short enough time scales, the defect occupation is frozen, and 
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FIG. 4. (a) Resonance fluorescence (0.25 nW at A = 951.6040 
nm, P — 0, 0. 1 s integration time) from a dot in Sample C with 
dcap = 150 nm (black points; red line Lorentzian fit), (b) Monte 
Carlo simulation using N 2 d = 1.0 x 10 10 cm~ 2 , T = 0.8 peV, 
N = 10, p = 4.4% (black points; red line Lorentzian fit). 



the dot's optical linewidth is close to transform-limited. The 
behavior for longer integration times, Fig. 3(c), is reproduced 
in the simulations with the same parameters but by increasing 
N, the number of hole configurations, in accordance with the 
integration time in the experiment, Fig. 3(d). 

A key conclusion for sample B is that local fluctuations of 
hole charges are responsible for the spectral fluctuations and 
an increase in the optical linewidths in time-integrated spec- 
tra above the transform-limit. The first step in reducing the 
optical linewidths is to take control of these holes. The hole 
density (Nh = pN^d) estimated in sample B at P = is 
very small, and our experiments have shown that Nh is not 
related to the weak resonant excitation. It is roughly consis- 
tent however with the p-type background doping of ~ 10 14 
cm~ 3 . Eliminating these holes completely may be challeng- 
ing. However, based on this new understanding, we have pur- 
sued the idea of reducing their influence by forcing the holes 
to adopt a position further away from the quantum dots. This 
reduces a hole-induced electric field at the location of the dots 
by ensuring a closer match between the electric field from an 
occupied defect and its image charge. Fig. 4 shows resonance 
fluorescence from a dot in sample C with an increased cap- 
ping layer thickness, d cap = 150 nm. There are two striking 
features. First, the linewidth has reduced to 1 .4 peW. The av- 
erage linewidth on sample C is 1 .60 peW with standard devia- 
tion 0.22 /ieV. Second, the fluctuations in Fig. 3(a) disappear. 
We attempt to reproduce this behavior in the simulations by 
keeping N2D, T, TV and Nh exactly the same as for sample 
B, changing only the capping layer thickness, Fig. 4(b). This 
results in a close-to-Lorentzian line with FWHM 1.1 peV, in 
very close agreement with the experiment. Quantitative un- 
derstanding of the valence band localization centers has there- 
fore been achieved. 
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